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Solid state quantum computing proposals rely on adiabatic operations of the exchange gate among 
localized spins in nanostructures. We study corrections to the Heisenberg interaction between lateral 
semiconductor quantum dots in an external magnetic field. Using exact diagonalization we obtain 
the regime of validity of the adiabatic approximation. We also find qualitative corrections to the 
Heisenberg model at high magnetic fields and in looped arrays of spins. Looped geometries of 
localized spins generate flux dependent, multi-spin terms which go beyond the basic Heisenberg 
model. 
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I. INTRODUCTION 



Scalable quantum computation proposals, particularly 
in solid state architectures compatible with existing mi- 
croelectronic technology, are of great potential impor- 
tance. The exchange gateA, based on the Heisenberg in- 
teraction between localized spins is a key concept under- 
lying several proposed quantum computer architectures 
in semiconductor nanostructures, where the spin of a lo- 
calized electron serves as the single qubiti^i&A In most 
spin-based quantum computer proposals single-qubit op- 
erations involve rotations of individual spins by external 
magnetic field pulses (essentially a modified electron spin 
resonance technique implemented on individual spins) or 
by some other techniques (e.g. local g- factor manipula- 
tion by external electric field pulses). The two qubit op- 
eration in solid state, spin architectures is projected to 
be achieved by the exchange gate operation. The ability 
of a tunable exchange gate (which enables the exchange 
coupling to change from zero to a finite value within in 
a "reasonable" gating time) in carrying out the "SWAP" 
operation leads to the implementation of the universal 
2-qubit CNOT gate which, along with single qubit gates, 
provide a universal set of quantum gates. The perceived 
advantages of solid state quantum computation are its 
scalability (because of the existing semiconductor micro- 
electronics infrastructure), the low decoherence rate for 
spin states (as compared for example, with charge states 
in semiconductors), and the feasibility of precise control 
over the local inter-electron Heisenberg coupling through 
the exchange gate. Among the disadvantages is the in- 
evitable presence of some spin decoherence due to the 
background solid state environment^ and difficulties in 
the measurement of single electron spin states. How- 
ever, recent experimental advances 6 have demonstrated 
single spin measurements of localized electrons embedded 
in semiconductors. 

In this work we provide a detailed and quantitatively 
accurate theoretical study of two complementary as- 
pects of solid state spin quantum computer architec- 
tures. First, we develop a theory for obtaining an ac- 



curate map of the low lying Hilbert space of the quan- 
tum dot based spin quantum computer architecture (the 
Loss-DiVincenzc- architecture in coupled GaAs quantum 
dots) in the presence of an external magnetic field. Our 
calculation of the Hilbert space structure of the coupled 
quantum dot (with one electron on each dot) system dif- 
fers from earlier work on the problem in the sense that 
it is essentially exact for our modelA^ Earlier work ob- 
taining the Hilbert space structure of the coupled dou- 
ble dot quantum computer system used perturbative ap- 
proximations akin to Hcitler-London or Hund-Mullikcn 
theoriesAiS As explained later in this paper we manage 
to carry out an exact numerical diagonalization of the 
two-electron coupled-dot interaction Hamiltonian by bor- 
rowing theoretical techniques which have been extremely 
successful in elucidating the ground and excited state 
properties of the strongly correlated fractional quantum 
Hall system. Second, we describe and discuss (again us- 
ing the exact diagonalization technique as the theoret- 
ical tool) a subtle topological feature of a certain class 
of solid state spin ("cluster") qubits, which have at- 
tracted considerable recent attention, where clusters of 
spins are cleverly aligned to serve as qubits (rather than 
just single spins) offering certain advantages in quan- 



tum computation 
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We show that in looped ge- 



ometries these spin cluster qubits have higher-oder spin 
interaction terms (3-spin, 4-spin, etc. terms) arising 
from chiral interaction terms which will have serious ad- 
verse consequences for quantum computation by provid- 
ing nontrivial flux dependence which must be taken into 
account: 15 Our finding of the chiral Hamiltonian in spin 
cluster qubits may have possible consequences for all spin 
quantum computing architectures (and not just the semi- 
conductor nanostructure based onesA&). Again, the exact 
diagonalization technique is the theoretical method we 
employ to demonstrate, validate, and quantify the chiral 
spin Hamiltonian. 

The common theme running through the two com- 
plementary topics studied in this work is a thorough 
investigation of the precise applicability of the Heisen- 
berg Hamiltonian as a description of the spin interaction 
in exchange gate quantum computer architectures. In 
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particular, we are interested in figuring out the limita- 
tions, and the constraints on the Heisenberg interaction 
model as the underlying Hamiltonian for exchange gate, 
two qubit entanglement in the presence of an external 
magnetic field. In this context we also want to know 
what, if any, (and how large) the correction terms are to 
the Heisenberg model description of the exchange gate 
quantum dynamics in the solid state spin-based quan- 
tum computer. This study is motivated by the fact that 
precise knowledge of the exact Hamiltonian controlling 
the qubit dynamics is absolutely essential in develop- 
ing quantum computing algorithms (and architectures) 
since in quantum computation the Hamiltonian itself de- 
termines the programming codes. We must, therefore, 
know the Hamiltonian precisely. 

To study modifications to the Heisenberg Hamiltonian 
in tunnel coupled single-electron quantum dots in an ex- 
ternal magnetic field we consider a realistic Hamiltonian 
involving an equal number of parabolic quantum dots 
and electrons interacting through the Coulomb interac- 
tion. We find deviations from simple Heisenberg behav- 
ior where expected, when the quantum dots are strongly 
coupled and under intense magnetic fields. We find, using 
a variational ansataii and exact diagonalization, that at 
high magnetic fields the two-electron orbital states may 
be characterized by a vorticity. We define vorticity as the 
number of zeroes attached to each electron in the many- 
body wave function through the term: [z% — Zj) p , where 
z = x + iy is a complex coordinate in the x — y plane 
and the integer p is the vorticity of the state. We find 
level repulsion between states with either even or odd 
vorticity. We identify the parameter regime required to 
maximize level repulsion among wanted and unwanted 
states ensuring adiabatic operation of the exchange gate. 

Using perturbation theory applied to the extended 
Hubbard model we find other qualitative modifications to 
the Heisenberg interaction among weakly tunnel- coupled 
quantum dots. We find that in spin clusters formed from 
loops of three or more spins, many-spin terms couple to 
external sources of flux. We focus, numerically, on the 
triangular configuration in particular. Here, we show 
quantitatively that a three spin chiral term couples to 
flux through the triangle^ We investigate both the low 
and high magnetic field regimes and model the effective 
spin Hamiltonian in terms of two parameters: The flux 
passing through the triangle and the vorticity (or effec- 
tive flux) attached to each electron. We go on to study 
four spin configurations where three and four body terms 
modify the Heisenberg interaction>iS We find that four 
body terms depend on the flux through closed loops as 
well. 

Implementations of exchange based-only quantum 
computation with the least overhead, from a quantum 
computing perspective, involve several spins interacting 
simultaneously through the Heisenberg interaction. Al- 
though, this is not necessary. One may consider al- 
gorithms involving no more than two simultaneously 
coupled spins or geometries which exclude closed loops 



of spins. In what follows we analyze a special case: 
looped geometries of simultaneously interacting single- 
spin quantum dots which, as we will show, necessarily 
involve flux dependence. 

In Sec. II we present a model of several lateral, single- 
electron quantum dots in an external magnetic field. The 
model establishes four regimes defined by parameters re- 
lated to confinement, the external magnetic field, and 
inter-dot separation. In Sec. IIA we relate the coupled 
quantum dot model to a fourth order spin Hamiltonian 
based on well known perturbation theories of the ex- 
tended Hubbard model. Here, we broaden this treatment 
to include magnetic field effects. We find many-body 
spin terms which couple to external magnetic fields. In 
Sec. IIB we describe a variational theory of the many- 
dot problem which accesses non-perturbative regimes of 
the parameter space. In Sec. Ill we discuss implications 
of the perturbative and variational treatments to qubit 
proposals involving two, three, and four coupled spins. 
In Sec. IVA we use exact diagonalization of a physically 
plausible, Coulomb Hamiltonian to explore the low en- 
ergy Hilbert space of two, coupled, single electron quan- 
tum dots. We find that the variational ansatz agrees 
with exact results. The results suggest that electrons in 
coupled quantum dots capture vortices of the many-body 
wave function to screen the Coulomb interaction, at high 
fields. In Sec. IVB we compare exact diagonalization re- 
sults of three, coupled, single-electron quantum dots with 
the spin Hamiltonian derived using perturbation theory. 
We show, numerically that chiral, three spin terms cou- 
ple the spin states of a decoherence-free subsystem to 
external sources of enclosed flux. We parameterize the 
magnetic field behavior of the lowest spin states with 
a spin Hamiltonian which depends on the enclosed flux 
and the number of vortices (or effective flu:»i£) attached 
to each electron. We conclude in Sec. V. 



II. MODEL 

We study the low energy Hilbert space of N, lateral 
quantum dots containing N electrons lying in the x — y 
plane with the following Hamiltonian: 
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We focus on GaAs. In which case we have an effective 
mass m* = 0.067m e , dielectric constant e — 12.4, and 
g-factor g* = —0.44. We work in the symmetric gauge 
with magnetic field B = Bz. S is the total spin. The 
single particle potential confines the electrons to lie in 
parabolic wells centered at TV positions, R;: 
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FIG. 1: Schematic showing four regions characterizing spin 
and orbital degrees of freedom in several coupled quantum 
dots defined by the parameters in Eq. Q at fixed confine- 
ment, ujo- The vertical axis depends on the inter-dot spacing, 
7?, while the horizontal axis depends on the ratio between the 
cyclotron frequency, ui c , and ujq. The region of interest for 
quantum computing (spin Hamiltonian) yields a spin Hamil- 
tonian dominated by Heisenberg exchange coupling, JijSi-Sj. 
Here, the higher orbital energy levels of the quantum dot 
have much higher energy than the exchange splitting, A> J. 
Above this region, R/a S> 1, the electrons in each dot do not 
interact strongly. At high magnetic fields and near R/a ~ I, 
the electrons capture vortices of the many-body wave func- 
tion to form mixtures with A > J. Below the dotted line the 
small separation between dots allows single dot behavior, and 
therefore level crossing, A < J. 



where u>q is a parabolic confinement parameter which 
may be compared to the cyclotron frequency, lu c — 
eB/m*c. The confining potential localizes the electrons 
at N sites with inter-site separation R — |Rj — R,-|. 
R may be compared to the modified magnetic length 
a = ^Jhc/eB(l + 4:UJq/uJc)~4 . Solutions of the above 
Hamiltonian take the form -4[|i/>) (g> \X)n] where and 
| A) jv are the orbital and spin parts of the wave function, 
A is the antisymmetrization operator, and the subscript 
denotes the number of electrons and quantum dots. 

The parameters in H define several regimes relevant for 
quantum computing architectures utilizing similarly con- 
fined single-electron quantum dots in an external mag- 
netic field. Fig. ^ depicts four separate pieces of the pa- 
rameter space with the confinement parameter, uo, fixed. 
The solid line encloses an area in which the excited, or- 
bital states of the quantum dots have high energy, A, 
and the inter-dot coupling between two dots maps onto 
the Heisenberg Hamiltonian, J12S1 • S 2 , as originally en- 
visaged in Ref. [l|. We show, using both perturbation 
theory in the Hubbard limit (Sec. IIA) and exact diago- 
nalization (Sec. IVB), that in this regime several coupled 
quantum dots, N > 2, involve symmetry breaking many- 
spin terms with nontrivial flux field dependence. Above 
the spin Hamiltonian regime (R/a 1) weak inter-dot 
tunneling yields a small exchange interaction. Below the 
dashed line, higher orbital states of the quantum dot have 
a particularly low energy and therefore mix with the low 
energy states of the quantum dot, A < J. This regime is 



characterized by magnetic field dependent level crossings 
among potential qubit spin states and unwanted higher 
energy levels of the double dot system. Between the level 
crossing and spin Hamiltonian regimes we find, using a 
combination of exact diagonalization (Sec. IV) and vari- 
ational techniques (Sec. IIB), that electrons in the dots 
capture vortices of the TV-body wave function to screen 
the strong Coulomb interaction. Here we find that the 
first excited state of the many dot system mixes vortic- 
ity leaving the unwanted excited states of the quantum 
dot somewhat higher in energy than the spin splitting 
between the lowest states, A > J. In Sec. II we discuss 
how these regimes pertain to qubit proposals utilizing 
coupled quantum dots. 



A. Perturbative Expansion 

We first seek an approximation to Eq. Q that quali- 
tatively captures the structure of the low energy Hilbert 
space in terms of on-site spin operators in the limit of 
weak inter-site coupling. Such a Hamiltonian may be 
used to define qubit gates among spins localized in neigh- 
boring quantum dots. We work in the single band, tight 
binding limit, a good approximation in the limit A> J. 
We also take the on-site Coulomb interaction to be much 
larger than the tunneling energy. As a result we may 
also restrict our attention to singly occupied states. As 
a first approximation we then obtain the extended Hub- 
bard Hamiltonian: 

+ V ^ "»,qtti+l,a / + g*MB B - ^ S i> ( 3 ) 

i,(a,a')SU « 

where cj Q creates a fermion at the site i with spin a and 
rii a = c\ a c ia . Hh incorporates the on-site spin opera- 
tors, Si = \c\ a (T aa ic ia , , where tr are the Pauli matrices. 
In the presence of an external magnetic field the tunnel- 
ing coefficients are complex: tjj = \tij | exp(27ri<I > i : j/<I > o)j 
where $0 = hc/e is the flux quantum. The magnetic 
vector potential generates the Peierls phase: 

*y = f A • dr, (4) 

where the integral runs along a path connecting the 
sites i and j. Working in the limit \Uj\/U -€1 1 and 
» e < 1 we can self-consistently confine our at- 

tention to the single occupancy states of the full Hilbert 
space. For large v the extended Hubbard term favors 
double occupancy and our approximation breaks down. 
We consider this regime variationally in the next section. 

We expand Hn by applying a unitary transformation 
exp(iK)HH exp(— iK), where K is an operator changing 
the number of doubly occupied statesii2ii2i22i2ii22 We ob- 
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tain, up to constant terms, the following spin Hamilto- 
nian: 

i i,j 

+ £ CiJAitC-Si-Si + ix-Sfc.Si + i) 

+ (-S J --S fc + i)(-S i -S / + J) 

(-Si-S fc +j)(-S J -.S i + j)]+,?^, (5) 
where we define the following flux independent factors: 
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and the following flux dependent factors: 
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The second sum in Eq. JSJ) contains the usual Heisenberg 
term if we define Jij — 2\Uj\ 2 /U . The third and fourth 
sums modify the Heisenberg term and depend on the 
lattice vectors r connecting a site to its neighbor. The 
fifth sum is a three site sum over chiral^ terms around 
distinct, closed loops (A), denoted Xijk = <Bi,i,fcSj • Sj x 
Sj,. The sixth sum includes four sites around distinct, 
closed loops (□). The coefficients in the last two sums 
depend on the flux enclosed by three site loops, &ijk, 
and four site loops, &ijki- H c g applies to any half-filled, 
single-band, singly occupied lattice in the presence of an 
external magnetic field, and in the limits i/f « 1 and 
v < 1. 

The magnetic field dependence of the tunneling ma- 
trix elements may be calculated directly using a Wannier 
basis formed from dot-centered, Gaussian, single particle 
states. An explicit form for \ti2\ 2 /U may be found in 
Ref.© with a slightly different confinement than the one 
defined in Eq. @ • Ref. contains two results relevant for 
our discussion. One first observes that, in the Hubbard 




FIG. 2: Two possible looped configurations of lateral quan- 
tum dots containng single spins in a magnetic field. The tri- 
angular configuration contains three tunneling channels with 
one loop formed by the vertices 1231. The square configu- 
ration contains, in general, six tunneling channels with seven 
disctinct loops formed by traversing vertices 1231, 2342, 3413, 
4124, 12341, 12431, and 14231. 



approximation, \ti2\ 2 /U decreases exponentially with in- 
creasing magnetic field or inter-dot separation. Second, 
one finds J12 ~ \ti2\ 2 /U > 0, for all magnetic fields and 
inter-dot separations. 

Eq. © incorporates three and four body spin terms. 
The three body chiral term splits the energy between 
states involving third order, virtual tunneling processes 
along and counter to the applied vector potential. The 
phase, 27r$yfe/<&0; is the Aharonov-Bohm phase gener- 
ated by the virtual current moving around the flux en- 
closed by the three site loop. The chiral term vanishes on 
bipartite lattices as a result of particle hole symmetry 22 
ft plays a particularly active role in triangular lattices, 
Fig.0 

While the chiral term vanishes in the zero field limit, 
four body terms survive for N > 3. Four body terms 
have been discussed in relation to coupled quantum dot 
systems in a four site, tetrahedral geometry^ In the 
tetrahedral geometry there are three terms in the sum 
over four site loops (three distinct CPs) yielding a par- 
ticle exchange symmetric, four body contribution to the 
spin Hamiltonian: ~ J2ijkleai^i ' Sj)(Sfc • S/). We may, 
similarly, consider a two dimensional square geometry 
with equal tunnelling between all sites. Here we find 
three distinct, four-site loops, Fig. |21 These add to 
give the term discussed in Ref. 10. In the absence of 
diagonal tunneling only the exterior loop, 12341, sur- 
vives, yielding the symmetry breaking four body term: 
(Si x S a ) • (S 3 x S 4 ) + (Si • S 2 )(S 3 • S 4 ), which, through 
Cijki , shows flux dependence. 

As the inter-dot separation decreases we leave the 
single-band, single-occupancy limit and a real current 
(rather than a virtual one) splits the energy spectrum. 
The spin splitting in this limit depends on first order, 
direct exchange processes rather than the second order, 
superexchange processes discussed above. Our anaylsis 
based on the extended Hubbard model may miss qualita- 
tive features in the low energy Hilbert space by poorly es- 
timating the coefficients of different spin terms. Nonethe- 
less, symmetry requires that, in the weakly-coupled lat- 
tice limit Eq. (JSJ) includes all possible scalars formed from 
spin operators, up to and including four body terms. 
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While the magnitude of these coefficient may be esti- 
mated accruately by expanding the Hilbcrt space, we seek 
to study the qualtitative physics associated with changes 
in sign. To go beyond our single band lattice model we 
account for qualitative changes in sign by modelling or- 
bital effects with a variational ansatz. 



B. Variational States 

We now discuss a set of variational states which model 
the low energy orbital states of Eq. in an effort to go 
beyond the Hubbard limit discussed above. To obtain an 
accruate variational wavefunction we examine the form 
of the exact wavefunction in two limits: the upper left 
and lower right corners of Fig. ^ We then construct an 
ansatz which connects both regimes. We begin with the 
simplest system, N = 2. It is analytically soluble in two 
extreme regimes: Two well-separated one-electron "arti- 
ficial atoms" and a two-electron artificial atom in a high 
magnetic field. The first case is trivial and consists of two 
well separated quantum dots (akin to two well-separated 
one electron atoms not in a molecular state) with one 
electron in each dot, fi/o> 1. In this case we may ignore 
the Coulomb interaction. The non-interacting ground 
state consists of degenerate singlet and triplet states. 

In the second soluble limit (a two-electron artificial 
atom) two electrons lie in one parabolic dot in a strong 
magnetic field. In this case we take uj c /<-^o 3* 1 and 
R = 0. Correspondingly, the relative and z-component 
of angular momentum commute with the Hamiltonian. 
At large magnetic fields we may project Eq. onto the 
lowest Landau level (LLL), giving: 

oo 

H | Wo <c c ,iJ=0 = lL z + ^ VmPm, (6) 
m=0 

where we define 7 = |(a/w 2 + Auj 2 — lu c ) and L z is the 
total angular momentum in the z direction. The second 
term represents the LLL Coulomb interaction, projected 
onto eigenstates of relative angular momentum, to, via 
the projection operator, P m . The coefficients, V m , are 
the Haldane pscudopotentials 24 which, for the Coulomb 
interaction, decrease with increasing to, at large to. The 
unnormalized eigenstates of relative angular momentum 
are: 

\m) = ( Zl - z % T exp ( ~ N ^ N2 ) . (?) 

where z — x + iy. It can be shown directly that, because 
there is no center of mass motion, the above wave func- 
tions are also eigenstates of L z , with eigenvalue to. Thus 
the set of states \m) form an orthogonal set of eigenstates 
of Eq. ©, with eigenvalues E m = jm+V m . The relative 
angular momentum of the lowest energy state depends on 
the parameters in 7 and the form of the interaction. For 
the LLL Coulomb interaction, in the artificial zero field 



limit, the lowest energy state has to = 1. Increasing 
B lowers the confinement energy cost, 7TO ~ m/B, and 
raises the Coulomb cost, V m ~ yjB/(m+ 1), thereby 
raising to by one. The transition from one eigenstate to 
the next occurs when E m = E m+ \ which, for uj c /luq 3> 1, 
occurs at magnetic fields: 

where V m = V m /(e 2 /4irea) and C = 4neh 3 / 2 uj 2 m*e~ 7 / 2 . 
For huj Q = 3 meV we find C[T 3/2 ] - 1.2. The states, \m), 
are symmetric (antisymmetric) with respect to particle 
exchange if to is even (odd). The total wave function, 
A[\ip) (8 I A) 2] , must be antisymmetric. Therefore |A)2 is 
spin singlet (triplet) for to even (odd). Here, the index 
to may be interpreted as the number of zeros or vortices 
attached to each electron, allowing us to assign a vorticity 
to each spin state. 

We now construct variational states which reproduce 
the exact results discussed above and the low energy 
physics of the intermediate, physical parameter regime 
as well, R/a ~ 1 and uj c /luo ~ 1. The composite fermion 
theory^ of the fractional quantum Hall effect offers an 
accurate variational ansatz describing two dimensional 
electron systems at high magnetic fields. A composite 
fermion is the bound state of an electron and an even 
number of quantum mechanical vortices of the many- 
body wave function. The corresponding orbital wave 
function isM^: if) = J(j) 1 where is a weakly interacting 
fermion state, J a Jastrow factor, and ip the highly cor- 
related state of electrons. In isotropic, spinless systems 
J attaches an even number of vortices to the fermions 
in the antisymmetric state <f> yielding an antisymmet- 
ric electron wave function. In anisotropic systems with 
additional quantum numbers one may bind an even or 
odd number of vortices to each particle while preserv- 
ing the antisymmetry of the overall wave function^SiS 7 . 
Applying the composite fermion ansatz to the Hamilto- 
nian studied here we take <j> to be the non-interacting 
ground state of Eq. (JTJ . A more rigorous, but technically 
demanding, approach for large systems should use the 
weakly interacting ground state instead. We also take 
3 = Uj<k( z j - z k) p , where p = 0, 1, 2, .., giving: 

% = Y[(zi-z k y<f>. (9) 

j<k 

This is our initial, high field solution of Eq. at arbi- 
trary N. ip P reduces to Eq. Q at R = 0. For R/a > 1 
the fermions in the state <f> become localized on each dot 
leaving J constant ~ RP . In which case ip p reduces to 
the limit of two independent electrons. 

To improve the variational states in the low field 
regime, lu c /ujq < 1, we note that the Coulomb energy 
cost may be lowered by mixing with higher energy states 
of the quantum dot thereby increasing the average inter- 
electron separation. In the R = limit, rotational sym- 
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metry requires the addition of states with the same an- 
gular momentum. This leads to the following trial states: 

i>v = life " Zfe ) p i 1 + ^ tflt ) & ( 10 ) 
j<k 

where the variational parameter (3 controls the amount of 
mixing with higher energy levels of the quantum dot. The 
total raising operators W = b\ + ... + b^ N and at = a\ + ...+ 
ajy act on the Fock-Darwin basis states centered between 
the dots. The single particle raising operators are given 
by: b\ = (z*/2 - 28 Zj )/V2 and a] = i( Zj /2 - 2d z *)/V2. 
The above variational states include mixing with higher 
energy states of the same angular momentum because the 
operator b'a' does not change the angular momentum of 
a Fock-Darwin state. They will be tested in Section IV. 



not necessarily met if we change the parameters in Eq. 
slightly because J has exponential B and R dependence 
at large B and R, respectively^ In fact, experiments on 
coupled quantum dots, while pushing for shorter gate 
times (and hence larger exchange energies), may indeed 
leave the border defined by the solid line in Fig. E 9 ' 28 
It is therefore important to understand the low energy 
Hilbert space of the coupled dot system when Jy > A. 
We will show, for N — 2, that the variational states 
discussed in Sec. IIB capture the magnetic field depen- 
dence of A. At large fields the variational states describe 
a bound state between electrons and vortices of the N- 
body wave function. We find the smallest A, A min , to 
occur when the vorticity (the number of vortices attached 
to each electron) of the first excited states mix to form 
an anti-crossing. 



III. SPIN BASED QUANTUM DOT-QUANTUM 

BITS 

Gate operations on single and multi-spin qubits rely 
on the adiabatic evolution of the spin state under the 
unitary time evolution operator defined in terms of the 
appropriate spin Hamiltonian (e.g. Eq. J5J): 

Texpf-tjf H eS (t)dtj \\(t' = 0))jv, (11) 

where T indicates time ordering and T the duration of 
a gate pulse. The qualitative spin physics captured by 
the spin Hamiltonian, Eq. J5J), therefore plays a crucial 
role in defining gates formed from coupled quantum dots. 
Concurrently, the orbital states, Eq. (JTJJJ, can be used 
to calculate the parameters in Eq. JSJ and their regime 
of applicability. In forming quantum gates out of cou- 
pled, single spin quantum dots we study modifications 
to the Heisenberg paradigm and its implications to qubit 
proposals in three different systems: A) Two strongly 
coupled quantum dots. B) Three simultaneously, weakly 
coupled quantum dots. C) Four simultaneously, weakly 
coupled quantum dots. 

A. Two Quantum Dots 

Solid state qubit proposals often make use of the 
Heisenberg exchange interaction between spins in neigh- 
boring quantum dots. The exchange interaction offers 
the potential for a universal set of quantum gates through 
the adiabatic operation of the exchange gate with* or 
without 14 single spin operations. Application of the ex- 
change gate to the two dot system will be an adiabatic 
process if the energy between the lowest, unwanted ex- 
cited state of the double quantum dot and the highest 
spin state storing quantum information is much larger 
than the exchange splitting, A 3> Jij ■ This condition is 
satisfied in the spin Hamiltonian regime in Fig. ^ but is 



B. Three Quantum Dots 

An accurate characterization of the double dot system 
allows us to define the appropriate parameter regime in 
which to study several Heisenberg coupled quantum dots 
and associated magnetic field effects. Two spin states 
of a three quantum dot structure can serve as an en- 
coded qubit. We first construct encodings which protect 
quantum information stored in many-body spin states. 
These encodings assume a noise operator which, as a 
demonstration, we choose to be collective or Zeeman- 
like. We then search for degeneracies in the set of states 
generated by these noise operators. The S — 1/2 sec- 
tor of the N = 3 system provides a simple example of a 
quantum dot, dccohercncc- free subsystem 29-30 . Consider 
three electrons confined to three quantum dots whose 
centers lie at the vertices of an equilateral triangle as 
shown schematically in Fig. [3 In this case a decoherence- 
free subsystem makes use of a fourfold degeneracy at 
B = to protect quantum information stored in the 
qubit defined by |A)3, where A = or 1, from fluctua- 
tions in the Zeeman energy. The four states are 3 ^: 

|A> 3 ® |-l/2) 

= ^(IUT>+^ A+1 im> + w 2 - A |TU>) 

|A> 3 <8> 1 + 1/2) 

= -L (| tn) +^ A+1 | TIT) +^ 2 ~ A | ITT)) ,(12) 

where to = exp(27ri/3). The second term in the tensor 
product denotes the total z-component of spin. 

Up to second order, Eq. JSJ) allows an encoding against 
Zeeman-like or collective noise. By collective noise we 
mean an interaction between spins and the environment 
which acts the same on all spins. By construction, the 
Zeeman term may alter the energy difference between the 
states in Eq. (|12|l with different S z , but not A. However, 
the chiral term in Eq. (|SJ) acts non- collectively. The chi- 
ral and Zeeman terms remove all degeneracies required 
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to construct a qubit immune to fluctuations in the per- 
pendicular magnetic field. Explicitly: 

Xi23|A> 3 = ^(2A-l)>/3|A> 3j (13) 

where #123 = (12tJ/L/)sin (27r$i 2 3/$o) m the case Jy = 
J and |iy | = t for all i and j. 

Following Ref.|32we may now, using Eq. (|5jl, construct 
a projected spin Hamiltonian which acts on the encoded 
basis states |A)3: 

Hn=3 = F 3 (AW,AW,AW) ■ S + ^B 1W S„ (14) 

where S is a pseudospin operator defined by projection 
onto two encoded basis states, \X)n hi our case. Fjv is 
a basis dependent, effective magnetic field which may be 
tuned through suitable manipulation of fcy and depends 
only on the coefficients of the two body terms in Eq. J5J . 
Fat may be calculated from these two body terms using 
the exchange operator: = (4Si • Sj + Iij)/2, where I 
is the identity operator. As apparent from Eq. (|14f) . X123 
yields an effective Zeeman splitting between the encoded 
basis states of the three spin qubit. In Sec. IVB we ver- 
ify numerically that the chiral term is actually sizable 
in the spin Hamiltonian regime of Fig. ^ We therefore 
arrive at a revealing inconsistency in seeking a decoher- 
ence free subsystem from a looped, three spin system. 
Part of our motivation for simultaneously coupling three 
spins was to remove the Zeeman term as a potential noise 
source. However, we have only enhanced the system's de- 
pendence on the external magnetic field by coupling the 
three spins in a loop. 

In the event that we may control the flux through the 
three spin system, the chiral term offers an additional 
one qubit gate. This term yields two advantages. The 
first stems from a comparison with single spin operations 
using localized magnetic fields. The three spin object en- 
compasses a larger area than a single spin and therefore 
eases constraints on locally applied magnetic fields used 
in applying single spin gates. Second, exchange-only en- 
coded universality schemes require a large overhead and 
extremely accurate application of the exchange gate to 
implement elementary algorithms. The chiral term may 
offer some relief form these constraints using algorithms 
which include the new, encoded Pauli-Z gate in Eq. (jTlfl . 



C. Four Quantum Dots 

We now turn to the case of four coupled spins, the 
lowest number of physical spins supporting a decoher- 
ence free subspacei 33 ' 34 We begin with four quantum dots 
containing four electrons coupled with equal tunneling 
\t{j\ — t, including diagonal terms. Fig. [21 shows a two- 
dimensional schematic. With equal tunneling we find a 
decoherence free subspace among two S — states cor- 



responding to A = and 1: 

ia> 4 = 1 mi) 

+ I UH)+^ A+1 l UU)+^ A+1 | iUT) 

+ cj 2 - x \ mi>+^ 2 - A i tut>. (15) 

Including all single and two-body spin terms in Eq. JSJl, 
these states show no explicit magnetic field dependence 
(excluding the magnetic field dependence of \Uj\ dis- 
cussed in Sec. IIA). As for the three-body term, the spin 
Hamiltonian must respect the inter-site exchange sym- 
metry inherent in the lattice. In the basis |A)4 we find: 

ijfceA 

= X (2A ~ 1} ^ ^ e ^|A) 4 , (16) 

where e^u is the four component Levi-Civita symbol and 
the sum excludes I = i,j or k. As expected, the sum 
vanishes with tunneling \tj,j \ = t for all i and j even in a 
uniform, external magnetic field. 

Four spin terms have a simple representation in the 
|A)4 basis. The last sum in Eq. in the case \tij\ = t, 
involves three sums over four site loops. Writing the four 
spin terms with the exchange operator, we find that they 
act as the identity operator in the basis defined by |A)4. 
In this case we have a simple, projected Hamiltonian: 

H N=4 = F 4 (A (1 \A (2 \A (3y )-S. (17) 

It is important to note that F4 depends only on coef- 
ficients from two body spin terms of the form Si • Sj. 
The three and four body spin terms (and therefore ex- 
ternal sources of flux through closed loops) do not affect 
the symmetric four dot structure with diagonal tunnel- 
ing. They must maintain inter-site exchange symmetry 
imposed by the lattice, precisely the symmetry exploited 
in constructing the decoherence free subspace, | A) 4 . 

We now consider symmetry breaking effects. In the 
absence of diagonal tunneling (ii.i+2 = 0) only the ex- 
ternal loop, around vertices 12341 in Fig. |3 in the last 
sum of Eq. J5J survives. The external loop alone breaks 
particle exchange symmetry. The additional term can be 
written: (Si x S 2 ) • (S 3 x S 4 ) + (Si • S 2 )(S 3 • S 4 ), ex- 
cluding two body spin terms. The entire looped term, 
including two body spin terms, contributes the following 
term to Hn =4 : — C\2?,a{S x + y/3S y ). From this term we 
see that in the square geometry fourth order terms not 
only modify the Heisenberg interaction, and hence the 
effective magnetic field F4, but also add and effective in- 
plane field in the |A)4 basis. The size of this effective, 
in-plane field depends on the real, external flux piercing 
the square plaquette through C1234 ~ cos(27r$i234/ < I > o)- 

Additional symmetry breaking occurs during gate 
pulses crucial to encoded universality schemes. In or- 
der to implement Pauli gating sequences on the encoded 
four spin qubit we must tune F 4 and therefore the tunnel- 
ing matrix elements iy . When applied to a decoherence 
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Single Qubit 




Encoded Hamiltonian 

F 2 -S 
F 3 -S + ^B- m S z 



F 4 • S — C1234 [S x + Sy^/a\ 



F 4 • S 



Spin Basis 

|A) a 
|A) 3 
A) 4 

|A) 4 



F 4 • S + f E B ljk e tjkl S z + d (t-yU s ) |A) 4 



FIG. 3: Table showing possible multi-spin qubits in the left 
column, the encoded Hamiltonian in the center column and 
the corresponding two-state basis in the right column. In the 
left column arrows denote single spins while lines indicate tun- 
neling channels which are all equal except in the last row. The 
Hamiltonians in the center column are written in terms of the 
encoded spin S defined via the two encoded basis states | A)jv- 
For N — 2 the basis states are the S z = 0, singlet and triplet 
states while for N = 3 and 4, |A)jv is defined in Eqs. 112H and 
IU5L respectively. The effective magnetic fields Fjv depend 
only on the coefficients of two-spin Heisenberg terms, A^\ in 
Eq. JSJ. The factors B and C depend on the coefficients of 
three and four spin terms in Eq. JSJ and therefore the flux 
through closed loops. 



generate an orthogonal set of Wannier basis statesASS 
The limited number of Wannier basis states allows for 
high accuracy only in a regime where the Coulomb in- 
teraction may be treated perturbatively. However, in 
our treatment we are able to access the strongly corre- 
lated regime by including up to ~ 10 5 Fock-Darwin basis 
states with z-component of angular momentum less than 
twelve. We use a modified Lanczos routine to obtain 
the ground and excited states. This technique yields the 
entire spectrum. However, here we focus on the low- 
est energy states. The energies for N = 2 converge to 
within 1/ieV upon inclusion of more basis states and may 
therefore be considered exact. While, for the N — 3 sys- 
tem, the ground and excited state energies converge to 
within 6/ieV upon inclusion of more basis states giving 
a strict variational bound to the accuracy. However, the 
slow convergence is due to corrections in the overall con- 
finement energy cost, ~ 10 meV. The energy differences 
quoted here converge much faster (< 2 ^eV) as we in- 
crease the number of basis states and may therefore be 
considered exact, with a few exceptions. These excep- 
tions occur near degeneracy points where our Lanczos 
routine requires a prohibitive number of steps to discern 
between two states whose energies are to within 5fieV 
of each other. In these rare, but important, cases we 
extrapolate between the nearest convergent energies. 



free subspace an encoded Pauli gate composed of Heisen- 
berg terms must, by construction, involve a spin spe- 
cific asymmetry. An example was considered in Ref. Il5t 
I £31 1 = |*23 1 = I £34 1 = t(l + S), where S is a number and 
all other \tij\ — t. The sum over chiral terms in Eq. i|lt)|) 
then induces an energy splitting 24irV3tJ5AB z /(U<i>o) 
between the states with A = and 1 for $yfe/$o <C 1. 
Here A is the area of the triangle denned by the vertices 
123 in the square geometry of Fig. This configuration 
is depicted in the last row of Fig. [31 The table summarizes 
five spin cluster qubit configurations and their encoded 
Hamiltonians written in the \X)n basis. From the table 
we see that configurations which break inter-site symme- 
try (rows two, three, and five) have non-Heisenberg terms 
which depend on the flux through closed loops. 



IV. NUMERICAL RESULTS AND DISCUSSION 

The accuracy of our perturbative and variational anal- 
yses may be checked numerically. We study two systems 
in particular, two electrons in two, adjacent quantum 
dots and three electrons in three adjacent quantum dots 
arranged in a triangle. We diagonalizc the full Hamil- 
tonian, Eq. QJ, in several regimes, including R/a ~ 1 
and uj c /uJo > 1. We construct the matrix representing 
H in the Fock-Darwin 35 * basis centered between the dots. 
Previous studies have employed diagonalization of simi- 
lar Hamiltonians using several dot centered basis states. 
This technique requires lengthy, numerical routines to 



A. Two Quantum Dots 

We seek a quantitatively accurate description of the 
boundaries and underlying physics of all regions depicted 
in Fig. 2] While we find that the perturbative expansion 
in Sec. IIA is valid for R/a > 1 and u> c /ujo < 3, the 
remaining portions of the parameter space involve long 
range correlations. Using the N — 2 system we check the 
accuracy of the variational ansatz discussed in Sec. IIB 
in several limits. We expect that the variational states 
discussed there should remain valid for N > 2, with ap- 
propriate modifications. 

We begin with the R = 0, lowest Landau level limit 
discussed at the beginning of Sec. IIB. The LLL approx- 
imation cleanly brings out the physics behind the high 
field spin transitions in two-electron quantum dots but, 
as we have discussed, needs modification at low magnetic 
fields. Fig. 0] plots E m — -E gn d versus B for the four low- 
est energy states, m — 1,2,3, and 4 with S z = 0. The 
parabolic confinement parameter is hug = 3 meV. Cusps 
appear at E m — -E gn d = where the ground state changes 
at B m signaling a change in the number of vortices per 
electron. (Note that the relation for B m , Eq. JSJl, is valid 
for lu c /luq S> 1.) The ground state clearly shows a num- 
ber of spin transitions with increasing magnetic field. 36 
Furthermore, the second highest excited state becomes 
degenerate with the third at level crossings which occur 
at magnetic fields between ground state transitions. This 
suggests that quantum information stored in the two low- 
est energy spin states in neighboring quantum dots be- 
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FIG. 4: Energy of the four lowest states of a single quantum 
dot (7? = 0) with two electrons under a perpendicular mag- 
netic field in the lowest Landau level plotted as a function of 
perpendicular magnetic field. The ground state energy is set 
to zero and the parabolic confinement parameter is huo = 3 
meV. The ground state, with orbital wave function given by 
Eq. f7p. alternates between spin singlet (S = 0) and triplet 
(S — 1) as a function of magnetic field. The spin singlet and 
triplet states correspond to even and odd angular momentum 
quantum numbers, m, respectively. 




B[T] 

FIG. 5: The top panel plots the energy of the four lowest 
states of Eq. Q obtained by exact diagonalization as a func- 
tion of magnetic field with the ground state energy set to 
zero. The separation between parabolic dots is now R = 10 
nm. The parabolic confinement parameters is huio = 3 meV. 
Transitions between spin singlet and triplet states remain. 
The bottom panel shows the overlap of the exact ground state 
and the trial states given by Eq. © . The number of vortices 
attached to each electron increases with magnetic field from 
p — to p = 3. As in Fig. 21 singlet (triplet) states correspond 
to even (odd) values of p. 



comes susceptible to leakage when the dots are brought 
very close together. 

We now turn to the case with finite inter-dot separa- 
tion, R > outside of the LLL. The top panel in Fig. [5] 
shows the four lowest energies obtained from exact diag- 
onalization of Eq. Q versus magnetic field. The energy 
zero is taken to be the ground state. We have chosen an 
inter-dot separation of R = 10 nm, confinement hojQ = 3 
meV, and S z = 0. The energy of the first excited state 
gives the effective exchange splitting which changes sign 
through successive spin transitions at each cusp. The 



results are qualitatively similar to the results shown in 
Fig. 0] but are entirely unexpected. Vortex attachment 
non-perturbatively lowers the Coulomb energy of uniform 
states but docs not necessarily apply to highly disordered 
systems. Yet, the intriguing oscillations in the effective 
exchange interaction seen in Fig. suggest just this and 
therefore require further study. 

In comparing Figs. 01 and [SJwe find further differences. 
At low fields, the top panel of Fig. [S] correctly shows a 
spin singlet ground state at B = rather than a triplet 
state as shown in the unphysical, LLL limit of Fig. 0] 
Most importantly, the degeneracies in excited states at 
5 = 0, 2.4, 5.2 and 8 T begin to lift, giving A min > J. As 
opposed to the level crossing in the single dot, R = case 
discussed earlier, the breaking of rotational symmetry 
forces an anti-crossing among the first and second excited 
states. At small to intermediated inter-dot separations, 
R/ a i; 1) the higher excited states are perturbed, single 
dot states with a nearly uniform charge density. 

A large anti-crossing among the two lowest excited 
states protects the quantum information stored in the 
entangled state of two strongly coupled quantum dots. 
Experimental uncertainties in R and cuo may eventually 
lead to the strongly coupled regime. Careful study of 
the states making up the anti-crossing is therefore cru- 
cial. The bottom panel of Fig. plots the overlap of the 
exact ground state and the variational state, Eq. ©, at 
R = 10 nm. Triplet (singlet) spin states correspond to 
odd (even) values of p, as in the R = case. The over- 
laps drop to zero when the particle exchange symmetry 
of the orbital wave function changes. We have checked 
by direct calculation of the density that, by B ~ 9 T, the 
modified magnetic length has become small enough to 
localize the electrons on each dot. The surprisingly high 
overlaps prove that vortex attachment is a valid ansatz 
even in the highly localized regime. At large dot sepa- 
rations, R > 40 nm, the Coulomb interaction lowers to 
a point where the splitting between spin states is near 
zero at large B. However, we have checked that even here 
the overlaps remain large. Another important feature of 
Eq. © is that the p = state does not take into account 
the Coulomb interaction. The overlaps near B = are 
correspondingly lower. 

The top panel in Fig. El plots the exact energy spec- 
trum, as in Fig. but for R — 20 nm. Here we see 
that, at large magnetic fields, the large separation be- 
tween electrons localized on each dot suppresses the ex- 
change splitting. However, several spin transitions still 
remain. The bottom panel in Fig. shows the energy 
of the four variational states, Eq. i|10|) . with p = 0, 1, 2, 
and 3. We take the ground state to be the zero in en- 
ergy. We obtain the energy by orthogonalizing the four 
variational states and diagonalizing Eq. £Q) in this four 
state basis. These variational states are an improvement 
over Eq. Jj|J). They include mixing with higher energy 
levels of the dots. The mixing is tuned with the vari- 
ational parameter /3. We minimize the energy with re- 
spect to (3 at each B. The parameter (3, of the ground 
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FIG. 6: The top panel shows the same as Fig. |5] but for a 
dot separation of 7? = 20 nm. The bottom panel shows the 
energy of the trial states in Eq. J I (It as a function of magnetic 
field. The ground state energy is set to zero. The energies 
are obtained by orthogonalizing the four modified variational 
states with p — 0, 1, 2, and 3 and diagonalizing Eq. Q. The 
variational parameter j3 is chosen to minimize the total energy. 




FIG. 7: Schematic diagram showing the area enclosing flux 
in the N — 2 and N = 3 systems defined by the density of 
the iV-body ground state. The dots represent the centers of 
the parabolic quantum dots. The patterned region for the 
N — 2 system defines an area A' symmetric about the axis 
joining the two dots. For the N = 3 system the triangular 
region encloses an area A. A correction to the triangular 
region, shown by the three patterned additions, defines an 
area « §A'. 



state varies from 0.02 at B = to 0.0006 at B = 5 T 
showing that large magnetic fields all but suppress Lan- 
dau level mixing. The exchange splitting obtained with 
the variational states compares well with the exact value. 
Furthermore, in the range B = 1 to 4 T, the second ex- 
cited state captures the essential features of the corre- 
sponding exact results. Rotational symmetry breaking 
forces the higher excited states to open an anti-crossing 
observed near B — 0, 2.4, and 4.3 T. The states at the 
anti-crossings in Fig. are similar to the states making 
up the level crossings in Fig.^J For example the electrons 
in the first excited state at 2.4 T in Fig. Inform a two and 
zero vortex mixed state in a 56% to 44% ratio, as opposed 
to the ground state which holds one vortex per electron, 
to within 98%. To evaluate the anti-crossing explicitly 
we note that for i?/a < 1 the asymmetry in confinement 
acts as a perturbation. We may rewrite the confinement 
potential up to an overall constant: 

V N=2 (r) = ^(\rf-\ X \R). (18) 

The second term breaks rotational symmetry and forces 
an anti-crossing among the lowest two excited states. It 
is important to note that the two lowest excited states 
involve states of even vorticity. Symmetry allows these 
two states to mix yielding an anti-crossing as one may 
find by diagonalizing the rotational symmetry breaking 
term in the even- vorticity subspace. The matrix elements 
are: m*uiQR/2 (ip p t\\xi\ + \x2\\ip p ), where, near B = 2.4 
T for example, p and p' may be or 2. These matrix 
elements give an anti-crossing A min ~ m*uj'^Ra + J . This 
is in contrast to ground state transitions between states 
with even and odd vorticity. Here the states ip p and ip p +\ 
cannot mix, allowing the exchange splitting to change 
sign. 



We stress that the top panel in Fig. is obtained by 
diagonalization of the Eq. with ~ 10 5 basis states 
while the lower panel is obtained by the same method 
but with four, physically relevant basis states. The agree- 
ment breaks down at larger fields, B ~ 5.6 T, because 
we have not included the p = 4 variational state in the 
excited states. Inclusion of variational states with large 
p is necessary at larger fields. The excellent agreement 
obtained thus far demonstrates that the plethora of spin 
transitions in strongly coupled double quantum dots orig- 
inates from a swapping of the particle exchange symme- 
try associated with vortex attachment. 

We may parameterize the high field effects of vortex at- 
tachment in an effective spin Hamiltonian based on the 
above numerical results and our analysis in Sec. IIB. 
Note that the exchange interaction changes sign in a 
roughly periodic fashion as each electron captures an ad- 
ditional vortex. The vortex may be interpreted, by equat- 
ing its Berry's^! phase to an Aharonov-Bohm phase, as 
additional flux^ The confinement, determined by loq, 
fixes the area defined by the electronic wave function, 
A', depicted schematically by the patterned region on the 
N = 2 side of Fig. [7| From Fig. \5\ we note that we may 
count the number of vortices attached to each electron us- 
ing the flux through the double dot system, BA' / $o- The 
parameter A 1 is fixed by requiring that the net Berry's 
phase swept out by one quasiparticle (the electron plus 
the attached vortices) circling the other quasiparticle as 
it encloses the double dot system, 2tt[BA'/^q — (N— l)p], 
must vanish for p = 2m& The net flux includes the effec- 
tive flux due to the Berry's phase associated with attach- 
ing p vortices to each electron. The data in Fig. [S] for 
example, show that at B — 5.4 T the flux through the 
double dot system exactly cancels the effective flux due 
to the attached vortices (two for each electron). By fixing 
A' in this way, nBA' /$o increases by integer multiples 
of 7r as each electron captures an additional vortex. We 
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may then write a parameterized spin Hamiltonian (up to 
second order in t 2 /U): 

H$ = 2J 12 S! • S 2 , (19) 

where J12 ~ ^-cos ( ^ . We determine A' from our 

numerical data and set t = \tn\- We have, for simplicity, 
excluded the Zeeman term. From Fig. [5] for example, we 
find the center of the p — 2 region to be Q>o/A' « 5.4T 
which gives A' ~ 800 nm 2 . The parameter A' suggests 
confinement of an appreciable part of the single electron 
density to within a radius of ~ 10 nm. When we insert 
the magnetic field dependence^ of \t\ and U into H^S 
we obtain qualitative agreement with our numerical es- 
timates of Jij at all magnetic fields. But, without the 
.2 

cosine term, jj remains positive for all B. 

B. Three Quantum Dots 

We now study the N — 3 system where the quantum 
dots lie at the vertices of an equilateral triangle with side 
lengths R = 40nm4^ We know from the previous section 
that a large inter-dot separation will prevent unwanted 
excited states of the quantum dot from approaching the 
spin states defining our qubit. We further expect the 
analysis of Sec. IIA to hold only for low magnetic fields 
while, at large fields, electrons capture vortices and ini- 
tiate spin transitions. As a consequence, an external 
magnetic field has three noticeable effects. 1). At low 
fields the length scale, a, is set by confinement and the 
flux enclosed by the triangular loop will dominate the 
magnetic field dependence of the states in Eq. 2). 
At higher fields the length scale shrinks with increasing 
magnetic field. The inter-dot tunneling matrix elements 
will be suppressed as the electrons become more localized 
on each dot. 3). The electrons will simultaneously cap- 
ture vortices to screen the increased Coulomb interaction. 
The latter effect, as for the N = 2 system, should, in the 
appropriate parameter regime, lead to oscillations in the 
total spin of the ground state as a function of magnetic 
field. 

Fig. [5] shows the energy of the lowest states obtained 
from exact diagonalization of H in the S z — 1/2 sector as 
a function of magnetic field. The confinement parameter 
is taken to be Hloq = 6 meV for the top panel and 3 meV 
for the bottom panel. The energy of the state |A = 0)3 
is set to zero. At B — the two lowest energy states 
have total spin S = 1/2 and are degenerate, as expected 
from the reflection symmetry of the triangular confining 
potential. The next highest state has S = 3/2 which 
corresponds to Q\tij\ 2 /U ~ 0.13 meV. Above this state 
we find (not shown) the higher excited states to lie above 
1 meV. 

We focus first on the low magnetic field data. As we 
increase the magnetic field the magnetic vector potential 
breaks the symmetry of the confining potential leading to 
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FIG. 8: Energy of the three lowest energy states versus per- 
pendicular magnetic field obtained from exact diagonalization 
of Eq. O in the S z = 1/2 sector for TV = 3. The centers of 
the lateral, parabolic quantum dots lie at the vertices of an 
equilateral triangle with 40 nm side lengths. The dotted line 
has total spin S = 3/2 while the dashed (|A = 1)3) and solid 
lines (|A = 0)3) have S = 1/2. The energy of the S = 1/2 
state corresponding to |0) 3 is set to zero. The top panel has 
a parabolic confinement parameter hcuo = 6 meV while the 
bottom panel has hcuo = 3 meV. The change in confinement 
changes the effective area of the system and, as a consequence, 
the position of degeneracy points between B = 1 T and 3 T. 



a splitting between the two lowest states. The splitting is 
linear in B, for small B, as in Eq. I|13l) . We expect such 
a simple behavior because the two body terms in Eq. (J5J 
have weak magnetic field dependence at low magnetic 
fields, through \tij\ 2 /U. Alternatively, the chiral term 
annihilates S — 3/2 states. Only two body terms in 
Eq. affect the S = 3/2 state. We therefore expect 
that the energy of the S — 3/2 state, Eg— 3/2, decreases 
very slowly with increasing magnetic field at low fields 
(while Eg= 3 /2 — £7|o) 3 should increase linearly). Here 
the contribution from the chiral term is sizable and is, 
for the parameters studied here, larger than the Zeeman 
splitting of a single spin in GaAs, w 0.025£?[T] meV. 

We now turn to the high field effects in Fig. |SJ For 
highly localized, non-interacting particles we expect the 
flux in £>i23 to be BA/$o, where A is simply the the area 
of the triangle defined by the centers of the three quan- 
tum dots. However, in our system, the parabolic confine- 
ment will not perfectly localize the interacting electrons. 
The area swept out by a closed loop around the bulk of 
the wave function will enclose an area larger than the 
triangle. Fig. [7| shows a schematic representation of the 
new, larger area encompassed by the N — 3 system. The 
additional area due to the expansion of the wave function 
is » 3A' /2, where A' is the area enclosed by an equivalent 
N = 2 system. We may use our analysis from the previ- 
ous section to determine the area added to the triangle. 
The net Berry's phase associated with virtual tunneling 
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processes of quasiparticles around the N — 3 system will 
be: 



„ BA WA' ,„ \ 

2tt 1 (N — l)p 

V$o 2$ 1 ; V 



(20) 



The additional flux will appear in the flux dependent 
factors in Eq. J3J). We rewrite the three spin Hamiltonian 
in the case of equal tunneling, t (excluding the Zeeman, 
extended Hubbard, and fourth order terms): 



24t 3 



U 2 



-sm 



2tt$ 



(3) 



Si • S 2 x S 3 



(21) 



where J ~ ^-cos (^pj. We define $( 2 ) = BA' and 

$ (3) s» B(A + 3A'/2). Note that integer p, in Eq. ffity. 
does not contribute to the chiral term. Therefore, vortex 
attachment does not directly affect the flux in the chiral 
term. Furthermore, the cosine in the first term param- 
eterizes large magnetic field behavior while the sine in 
the second, chiral term, was derived using perturbation 
theory. H cff ' allows us to predict the degeneracy point 
of the three spin term using the degeneracy point of the 
two spin term with the following energies: 



Es=3/2 - -B|0> 3 



^|1>3 _ E \0) 3 



6t 2 /7T$ (2 ) 
COS — 

U \ $0 



6\/3i 3 . 
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From the dotted line in the top panel of Fig. [8] and 
Eq. ([HJ) we find $ /4A' w 2.2T. Using Eq. |g3$, we pre- 
dict the dashed line to cross the x-axis near B 1.5T, 
where we take A to be the area of an equilateral triangle 
with 40 nm side lengths. A similar analysis yields good 
agreement for the bottom panel of Fig. [3] 

The slope of the energy splitting between the two low- 
est states in Fig. |S] allows us to estimate t/U for this 
system using Eq. I|22|) . t/U is largest and only weakly 
magnetic field dependent at low B. Taking A+ |A' from 
above and v = 0, we obtain t/U ~ 0.09 for the top panel 
and 0.19 for the bottom panel which shows that our ex- 
pansion in t/U is consistent. For N = 3, only odd pow- 
ers of tij allow linear magnetic field dependence in the 
splitting, showing that, excluding double occupancy, the 
magnetic field dependence captured by the chiral term is 
accurate upto!?^|ij. 



V. CONCLUSION 

We show in this work that the Heisenberg model de- 
scription of the qubit coupling in the quantum dot spin 



quantum computer architecture applies only in a limited 
regime of the parameter space. In the GaAs quantum dot 
exchange gate architecture, the Heisenberg spin Hamilto- 
nian description applies only in the intermediate regime 
R/a > 1 and lu c /luq < 1. Using the exact diagonalization 
technique one can map out the precise low lying Hilbcrt 
space, and consequently, use this information in the de- 
sign of the quantum computer architecture. We have also 
discussed an interesting and non-trivial level-crossing pe- 
riodicity in the singlet-triplet energetics. Precise knowl- 
edge of the associated, low-energy Hilbert space could, 
in principle, be used to protect quantum information en- 
coded in the electron spin. 

Generalizing our exact diagonalization technique to 
spin cluster qubits formed by a two-dimensional array of 
electron spins localized in tunnel- coupled quantum dots, 
we show that the chiral term associated with the quan- 
tum phase picked up by an electron enclosing the mag- 
netic flux through closed loops must be included in the 
spin Hamiltonian. The existence of the chiral term in 
the looped spin cluster qubits modifies the Heisenberg 
interaction and is in some sense a decoherence mecha- 
nism for the simple exchange gate architecture (since the 
two-qubit SWAP operation is no longer determined by 
just the Heisenberg exchange Hamiltonian). We show 
in this paper how precise knowledge of the cluster ge- 
ometry, combined with exact diagonalization, provides 
us with the multi-spin Hamiltonian which would be re- 
quired for quantum computation with two or three di- 
mensional spin cluster qubits. Strictly one dimensional 
spin cluster qubits, which do not have any topological 
looping, have a small (but non-zero) chiral contribution 
and are therefore described, for the most part, by the 
Heisenberg Hamiltonian in the appropriate subspace of 
magnetic field, confinement, and dot geometry parame- 
ters. 

We emphasize that in this article we have considered 
a relatively simple model, defined by Eqs. (JIJ and J2J, 
for determining the applicability of the Heisenberg inter- 
action in describing the exchange gate operation. Dif- 
ferences in the confinement potential may change some 
of the quantitative aspects of our results but as long as 
the confinement consists of smooth potential wells, there 
should be qualitative agreement. The key issues we have 
addressed in this work is the regime of validity of the 
Heisenberg exchange gate operation in coupled semicon- 
ductor quantum dot quantum computer architectures as 
appropriate, for example, in GaAs-based quantum dot 
systems. In practice, we have obtained the conditions 
and constraints necessary for a coupled qubit system to 
behave as a coherent molecule as opposed to two decou- 
pled atoms. Adiabatic tuning between these two regimes 
enables the swap idea underlying the exchange gate. Of 
course, the issue of He atom-to-molecule transition in 
coupled quantum dot systems as well as our discussion 
of the exchange oscillations in the coupled dot system as 
a function of the applied magnetic field have implications 
beyond quantum computation. For example, a direct ex- 
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perimental observation of the exchange oscillations is of 
interest in quantum dot physics. 

Finally we mention that there are many other factors 
beyond the scope of our work (i.e. beyond the model 
defined in Eqs. and J5J) which affect the operation 
of the exchange gate. We cite three such example of 



recent interest which have been considered in the litera- 
ture: Inhomogeneous magnetic field effects^, spin-orbit 
coupling 40 —, and multi-valley quantum interference 4 ^. 

We would like to thank J.K Jain and K. Park for many 
helpful discussions. This work is supported by ARO- 
ARDA and NSA-LPS. 
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